Fluid mixtures of parallel hard cubes 
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The direct correlation function of a fluid mixture of parallel hard cubes is obtained by using 
Rosenfeld's fundamental measure approximation. This approximation is thermodynamically con- 
sistent (compressibility and virial equations of state are equal) and predicts a spinodal instability 
of the binary mixture for large-to-small side ratio larger than roughly 10, in qualitative agreement 
with simulations on the lattice version of the model. In two dimensions the system never demix, 
also in agreement with the simulations. 

PACS: 61.20.Gy, 64.75.+g 

The theory of liquids is still looking for its 'Ising 
model', i.e. an exactly solvable model which provides de- 
tailed information on the mechanisms involved in liquid 
behavior, and against which approximate theories could 
be tested. Instead the situation is right the opposite, 
and the only analytical knowledge we have about liquid 
models has been obtained within the Percus-Yevick (PY) 
closure approximation of the Ornstein-Zernike integral 
equation, and consists of the direct correlation function 
(DCF) of a fluid mixture of hard spheres with an ar- 
bitrary number of component species 0, and that of a 
one-component fluid of adhesive spheres [0 . No other ap- 
proximation has provided any new analytical result yet. 
Recently a very promising approximate scheme has been 
proposed: the fundamental measure theory (FMT) @, 
which, as in the scaled particle theory, is based upon an 
interpolation between the low and high density limits of 
the fluid. By an appropriate choice of geometric mea- 
sures, the FMT provides the following expression for the 
DCF of a fluid mixture of hard convex particles: 

<V( r ) = [Xo + Xi-R^(r) + X2>V(r) + X3^ ll ,(r)]/ AIl/ (r), 



(1) 

where V^ v {v), S^t) and i? M „(r) are, respectively, the 
volume, surface, and mean radius of the overlap region of 
particles of species /x and v whose center-to-center vector 
is r; /^(r) is the Mayer funcion, and the \j s are density- 
dependent coefficients which can be obtained by means 
of the scaled particle prescription. When applied to a 
mixture of hard spheres, Eq. (^) reproduces the PY solu- 
tion; besides, the same scheme can be generalized to any 
odd dimension and thus obtain the DCF of mixtures of 
hyperspheres (in one dimension it reproduces the exact 
DCF). In spite of this remarkable results, the FMT has 
never been applied to any molecular shape other than 
spheres, and it has recently been shown that, in general, 
Eq. (Q) cannot be obtained for arbitrary molecular shapes 
without further approximations p|. 

While the one-component hard-sphere model has pro- 
vided an excellent starting point for the study of sim- 
ple liquids, either as a reference system for perturbation 



theories (S) or as an ingredient for density funcionals of 
inhomogeneous liquids || , the situation is far poorer for 
additive fluid mixtures, due to the fact that the PY DCF 
predicts stability of a binary mixture of asymmetric hard 
spheres for any diameter ratio and any densities of the 
components [jjj . This places the demixing transition out 
of the scope of an analytical treatment. Furthermore, 
this fact settled the belief that, contrary to what happens 
in other phase transitions (such as freezing), demixing 
cannot be driven by a pure entropic mechanism. This 
conclusion was questioned by Biben and Hansen , who 
found, by means of a numeric solution of the Ornstein- 
Zernike equation with a more accure closure relation, 
that a mixture of additive hard spheres of diameter ratio 
larger than 5 and similar concentrations demixes. Fur- 
ther evidence of this fact has subsequently been provided 
by means of alternative approximations || . The question 
of wether a purely entropic mechanism can account for 
phase separation was affirmatively answered by Frenkel 
and Louis l| by providing the exact solution of a two- 
dimensional lattice model of two types of particles. For 
continuous three-dimensional models the controversy is 
still open. Simulations cannot yet help because for too 
asymmetric mixtures relaxation times are prohibitively 
large 0. Nevertheless, a recent simulation of a three- 
dimensional lattice model of hard cubes |l(| provides 
strong evidence of phase separation for mixtures of cubes 
of side ratio 3, while no such an evidence exists when the 
side ratio is 2 (in both cases the side of the small cubes 
is 2 lattice spacings). Besides its computational simplic- 
ity, this model is interesting in its own because even in 
the mixed phase in can completely fill the whole space at 
close packing, and therefore there is no trivial volume- 
driven phase separation. For the continuous version of 
the model, the fourth order virial expansion of the pair 
correlation function diverges at contact when the side ra- 
tio goes to infinity ficfl , thus suggesting the existence of 
a demixing transition. 

In this letter I will show how the FMT provides an ex- 
pression for the DCF of a mixture of parallel hard cubes. 
By means of this function the free-energy of the mixture 
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can be computed, and its concavity requirement leads to 
a spinodal instability for large enough side ratio. 

The starting point of the FMT is the choice of the ge- 
ometric measures. This can be easily done if the Fourier 
transform of the Mayer function can be expressed as a 
sum of products of single particle functions. In our case, 
the Mayer function of two parallel hard cubes of sides 
and a v is simply 

/^(r) = -0(0^ - |z|)0(<v - |y|)e(tv - |z|), (2) 

0(u) being the Heaviside step function (= 1 for u > 
and = otherwise), and — (cj m + a v )/2. Therefore, 
since the Fourier transform of any of the step functions 
appearing in (0) is 



dx e 



ffj,(k) = 2fc _1 sin(cr M fc/2) and C M (fc) = 2 cos(cr Al fe/2), then 
the Fourier transform of the Mayer function can be ex- 
pressed as 

iUk) = 1 (4 3 )(k)4°)(k) +4°)(k)4 3 )(k) 

+w< 2 ) (k) • w« (k) + w« (k) • w( 2 ) (k)) , (3) 

where 

«>l 3) (k) =r p (fc x )f Al (fc y )f At (fc z ), (4a) 



w^ 2) ( k ) = [C^{kx)T^k y )T^{k z ),T^(k x )C^{k y )f^(k z ), 



Tu\kx}Tu{ky)C,U,{kz) 



(4b) 



w«(k) = [^(k x )Uk y )Uk z ),Uk x )Mk y )Uk z ), 



C,^Q i x)C,^{k y )T^(k z ) 



W^(k)^Uk X )Uky)Uk Z ). 



(4c) 
(4d) 



(the numerical coefficients have been chosen for conve- 
nience). The uniform limit of 7ij(r) is £j (i = 0,3), and 
that of nj(r) is Cj'u/3 (j = 1, 2), where u = (1, 1, 1) and 
the £s are the average density, mean radius, surface and 
volume of the cubes: 



(e ,a,6,e 3 ) = E( 1 'Y' 6(7 '' CT 



(7) 



The ith averaged density [i — 0, 1, 2, 3) is a function with 
dimensions (volume)' 1-3 )/ 3 ; therefore whose dimen- 
sions are (volume) -1 , can be expanded as 



$ = an + Anin 2 + Bn 2 n 2 n 2 , 



(8) 



where a, A and B are, respectively, a scalar, a 2-index 
tensor and a 3-index tensor, all dependent on 71,3. The 
scaled-particle prescription imposes on $ the following 
differential equation la]: 



4> 



<9n 3 



(9) 



which yields a = a,Q — ln(l — n 3 ), A = (1 — ?i 3 ) -1 Ao and 
B = (1 — n 3 ) -2 B , with a , A and B denoting a con- 
stant scalar, a constant 2-indcx tensor and a constant 
3-index tensor. The vanishing of the excess free-energy 
per particle in the zero density limit leads to so = 0. 
On the other hand, the uniform limit of this free-energy 
imposses a simplified form for the tensors: 



(Ao)ij = aiUjUf + 71%, 
(B )yfc = a^UjUjUfc + j2VLiSij, 



(10a) 
(10b) 



Now that we have identified the so-called fundamental 
measures, the FMT postulates the following excess (over 
the ideal one) free-energy functional of the set of density 
profiles for each component [p|: 



/3F CX [{p M (r)}]= J dr$(Mr),n,(r)}) 



(5) 



(f3 1 = fcsT) where <5 is assumed to be a function of the 
following weighted densities: 
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)iw«(r-r')dr', 


(6b) 
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)wf(r-r')dr' 


(6c) 
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)™f (r-r')dr', 


(6d) 



ai i2 and 7i j2 being numerical coefficients. Then 

ainin 2 + 71 ni • n 2 



-n ln(l - n 3 ) + 



Oiin 2 + 7 2 n 2 n 2 • n 2 



1 - n 3 



(l-n 3 ) 2 



(11) 



where rij = • u (i = 1, 2). 

Now we have to routes to obtain the DCF of the fluid 
mixture: (i) from Eqs. (g), (§) and (0), as |§ 



c M „(r - r') 



5F™ [{ Pm (t)}] 



(12) 



{p tl ( r )}={Pn} 



and (ii) from Eq. ([!]) , by explicitly obtaining the involved 
functions. This latter expression can easily be computed 
because the overlap region of two parallel cubes of sides 
(jfj, and er„ is a parallelepiped of sides L^ v {x), L^ v {y) 
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and L^ v {z) (provided they do overlap, i.e. /^i/(r) = — 1), 
where 



_ / 0>„ - Xfj, v if \s\ < 

a^ v - \s\ if |s| > X^y 



(13) 



(Xftu = |cr M - cr^ |/2). Therefore 



V^„(r) = L^ u (x)L fJ , l/ (y)L^ u (z), (14a) 
^(r) = 2{L lxv (x)L llv (y) + L flv (x)L tlu (z) 

+L^ v (y)L^ v {z)}, (14b) 

R^u(r) = - {i^(x) + L„ v {y) + L^{z)} . (14c) 

By matching both expressions of the DCF it follows that 
oti = 0, 7i = 3, a 2 = 5/432 and 72 = -1/48. Then 
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Xo = 
Xi = 

X2 = 
X3 = 
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-6) 2 



i CI 



266 i Cf 



(i-6) 2 (i-6) 3 36 (l-^) 4 



(15a) 
(15b) 
(15c) 
(15d) 



and the excess free-energy per unit volume turns out to 
be, in the uniform limit. 



$ = _£ ln(l-&) 



66 i 6 3 



216 (1-6) 2 ' 
The compressibility equation of state 
,<9P 



(16) 



(17) 



leads to 

f3P c 



6 66 i Cf 
1-6 (i-6) 2 108(1 -6) 3 ' K } 



On the other hand, the virial theorem provides a different 
expression for the pressure: 

PP V =& + \Y,P»P» j dvvVf^{v)g^{v), (19) 



with g^ v {v) the pair correlation function. Straightfor- 
ward manipulations lead to 

pp v = 6 - 4 PvPv G vv I d v I dz c ^(o-^ v , y> z )i 
»,u J ° J ° 

(20) 

from which it follows that P v — P c . The expression 
found for the DCF of this system has thus the remarkable 



property of being thermodynamically consistent without 
having impossed such a feauture anywhere in the ap- 
proximate scheme. Furthermore, this (unique) equation 
of state matches the second and third virial coefficients, 
which can be readily computed to be B 2 p 2 = 66 + §66 
and B 3 p 3 = 66 2 + 2666 +6 3 /108 (where p = YLP* = 
Co)) but this is implicitly assumed in the theory [pj. 

A direct comparison with the simulations of Rcf. |h]] 
is meaningless, because the latter are performed on a 
lattice system with a lattice spacing comparable to the 
small particle size; moreover, there are no other results 
of this kind available for this model, so the accuracy of 
the result just presented has still to be tested. For the 
one-component system, however, both the exact virial 
coefficients and the virial expansion of approximate the- 
ories (such as PY) are known up to seventh order |jTl| . 
From them we can see that the values obtained from Eq. 
( |l8| ) neither match the exact ones nor those obtained 
from the PY approximation, from the fourth coefficient 
on. Thus we can conclude that the present results are 
— as expected — approximate, but they differ from the 
PY approximation, contrary to what happens for hard 
spheres (the equivalence for hard spheres might indeed 
be an exception rather than a rule). 

But the most striking consequence which follows from 
the results presented thus far is the prediction of a spin- 
odal instability of a binary mixture of hard cubes for 
large enough side ratio. Furthermore, this spinodal can 
be analytically computed within this approximation. Let 
us see how it can be done. The instability is caused by 
a violation of the concavity of the free-energy. The con- 
cavity of this function can be expressed as the possitive 
defimteness of the hessian matrix of the free-energy per 
unit volume, / = F/V, as a function of the partial den- 
sities, p v . Since / is given by / = PvPv — P, with p v 
the chemical potential of the v component, by using the 
Gibbs-Duhem relation we can obtain [El: 



M -R ° 2 f R d P» 1 

M v x = P-R — 5 — = P-R — = — 
dp v dp\ dpx Pv 



5ux-c»x(0). (21) 



From the expression for the DCF the 2x2 matrix M 
of the binary mixture can be readily computed, and the 
stability condition follows from the possitiveness of its 
determinant — since all its elements are possitive. By 
defining r\ v = <J^p v , the packing fraction of species v; 
i] = r/i + r/2 = 6) t ne total packing fraction of the fluid; 
r = eri/<72, the large-to-small side ratio (r > 1); and 
x = 771/ 77, the relative packing fraction of the large com- 
ponent, a tedious but straightforward calculation yields 



P\P2\M\ 



>r 



1 



1 



3(r-l) s 



x(l — x) 
(22 



and thus the mixture is stable provided the expression 
within square braquects is possitive. 
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An immediate consequence of this result is that the 
mixture is stable if and only ifl < r < 5(1 + 
y/l - 1/25) ks 9.98. Above this threshold there is a 
region where the mixture phase separates (see figure), 
limited below by the curve 




0.0 1 1 1 1 1 

0.0 0.5 1.0 

X 

FIG. 1. Demixing phase diagram for a binary mixture of 
parallel hard cubes: rj is the total packing fraction and x the 
relative packing fraction of the large cubes. Curves represent 
the spinodal lines for side ratios 10, 11, 12, 15, 20 and 50 
(from top to bottom). 



which is logically missing. Again the equation of state 
reproduces the three first terms of the virial expansion, 
as expected from the FMT. The novelty is that it can be 
shown that the binary mixture of hard squares is always 
stable, whichever the size ratio. It is interesting to note 
that the authors of Ref. |l0| arrived to the same conclu- 
sion in their simulations. Full details on the hard-square 
mixture will be given elsewhere ]l2| ] . 

In summary, I have provided in this letter the DCF 
of a new model, namely a mixture of hard cubes, ob- 
tained by means of the FMT, an approximation intro- 
duced by Roscnfcld and thus far applied only to repro- 
duce the results of a mixture of hard spheres. Apart from 
the importance of having analytical results for new mod- 
els of liquids, the relevance of this calculations relies on 
the fact that a spinodal instability is predicted for the 
binary mixture, in qualitative agreement with the sim- 
ulations performed on an equivalent lattice model. The 
calculations may serve as a reference for further studies 
on related system, for instance, the orientational freezing 
of liquid crystals, if the model is extended to a mixture 
of parallelepipeds, or the polymer collapse, of which Ref. 
pPf reports some simulations. Of course it would also be 
interesting to apply this results to the lattice system, in 
order to see whether they can reproduce the simulations 
quantitatively. 



Thus we have achieved a qualitative agreement with 
the simulations of Ref. Jio| ; however, a direct comparison 
of the results is, perhaps, too naive. In the simulations 
the authors find demixing for r > 3, while this calcula- 
tions shift this value up to « 10. But the simulations 
are performed on a lattice, and then the size of the small 
particle (2 lattice spacings in this case) is also important, 
for a larger size with the same size ratio implies a larger 
number of accessible sites on the lattice for the small 
particles, with the corresponding gain in entropy. There- 
fore, if the small particles are larger, it is plausible that 
the stability of the mixture increases to larger size ratios. 
The results for the continuous model we are dealing with 
in this letter would correspond to the limit <j\ — > oo and 
(72 — > oo, while keeping o\j<Ji constant, what explains 
its larger stability. Some simulations could ascertain this 
conclusion. 

The two-dimensional case, i.e. a mixture of hard 
squares, is a simple reproduction of the calculations here 
presented. This is something unusual in view of what 
happens for hard spheres ||, for which the FMT can 
be applied only in odd dimensions. For cubes the cal- 
culations can be done in a similar way no matter the 
dimension. In two dimensions, for instance, the form of 
the DCF is that of Eq. p) except for the volume term, 
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